Heat radiation from long cylindrical objects 
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The heat radiated by objects small or comparable to the thermal wavelength can be very different 
from the classical blackbody radiation as described by the laws of Planck and Stefan-Boltzmann. 
We use methods based on scattering of electromagnetic waves to explore the dependence on size, 
shape, as well as material properties. In particular, we explore the radiation from a long cylinder 
at uniform temperature, discussing in detail the degree of polarization of the emitted radiation. If 
the radius of the cylinder is much smaller than the thermal wavelength, the radiation is polarized 
parallel to the cylindrical axis and becomes perpendicular when the radius is comparable to the 
thermal wavelength. For a cylinder of uniaxial material (a simple model for carbon nanontubes), we 
find that the influence of uniaxiality on the polarization is most pronounced if the radius is larger 
than a few microns, and quite small for submicron sizes typical for nanotubes. 
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I. INTRODUCTION 

Thermal radiation lies at the heart of modern statisti- 
cal physics and goes back to the beginnings of quantum 
mechanics more than a century ago pQ. Planck's law de- 
scribes the intensity T (per unit surface area and solid 
angle) of radiation of a black body at temperature T 



Here, u> is the angular frequency of radiation and h and c 
are reduced Planck's constant and the speed of light, re- 
spectively. Integration over angles and frequencies yields 
the well known Stefan-Boltzmann law[2] for the total 
power H radiated per unit area A, 

H/A = oT A , (2) 

wither = TT 2 k 4 B /{60h 3 c 2 ). 

Only recently, various phenomena, leading to modifi- 
cations of these laws have been explored. For example, 
the effect of spatial coherence of emitted heat radiation 
was studied by many authors [3HB] as this effect can be 
used in many technological applications [TUlOj . such as 
thermophotovoltaic and high-efficiency light sources. 

For real materials, Eq. (IT]) can be adjusted by intro- 
ducing the (angle dependent) emissivity of the material, 
which is unity for a black body. Considering objects with 
sizes smaller or comparable to the thermal wavelength 
At = hc/ksT (approximately 7.6 fim at room temper- 
ature), the radiated energy differs from the above equa- 
tions because of interference effects of the object with 
the emitted radiation. In other words, the emissivity 
then depends on the size and shape of the object. Ad- 
ditionally, if the object is smaller than the penetration 
(skin) depth, the emitted power is proportional to the 
object's volume, rather than surface area. The heat ra- 
diation of small spherical objects including these effects 
have been studied since the 1 9 70s [IT] [12] ■ Also, effects 
of excitations [13] and electric currents [14] on the radia- 
tion have been studied. Furthermore, recent studies on 



superscattering properties of subwavelength nanostruc- 
tures (e.g. nanorods) [T5] make such systems potential 
candidates for efficient heat transfer sources. 

Experimentally, the radiation of thin cylindrical ob- 
jects with thickness in the range of the thermal wave- 
length is very well accessible, and has e.g. been stud- 
ied using metal wires with interesting findings: 50 years 
ago, it was discovered that the radiation of a hot thin 
metal wire is significantly polarized [16j . Polarizations of 
28% [TB] and 50% [T7] in the direction perpendicular to the 
wire were measured for thin incandescent tungsten and 
silver wires respectively. In both studies the thickness of 
the wires was larger or comparable to the thermal wave- 
lengths. These findings triggered a number of studies on 
the properties of thermal radiation of sources of various 
designs, including platinum microwiresjlS, 1.9^, semicon- 
ductor layers in external magnetic fields [20], bundles of 
carbon nanotubes [2"Tl [2"2"] and SiC lamellar gratings [23 . 

For wires with thickness smaller or comparable to the 
thermal wavelength, the radiation was found (e.g. for 
platinum) to be polarized in the direction parallel to the 
wire becoming fully polarized as the width approaches 
zero [HI . Polarization effects have also been observed 
for radiation of bundles of carbon nanotubes [21] and are 
considered a simple way of finding the degree of align- 
ment inside the bundle. For carbon nanotubes, an ex- 
planation for this polarization, taking into account the 
electronic structure of the tubes, has been discussed in 
Ref. [22]. 

Recent work on the heat radiation of thin metal wires 
[24] provides experimental as well as theoretical results, 
albeit restricted to emission perpendicular to the cylin- 
drical axis. Also, a series of works [2"5T - f2"T] discuss ra- 
diation emitted by individual incandescent carbon nan- 
otubes. In Ref.[26 , the polarization of the radiation is 
studied both experimentally and theoretically using a 
model based on Mie theory. 

In this paper, we provide the general formalism to 
compute the heat radiation of arbitrary objects in terms 
of their classical scattering properties, which is part of 
the general framework [28] for non-equilibrium electro- 
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magnetic fluctuations involving multiple objects and ar- 
rays [29 35 . Thereby we give a more extensive deriva- 
tion of the corresponding results presented in Ref . [2"5] . 
which is also more general in terms of material proper- 
ties: We include the possibility of dielectric or magnetic 
losses, locality or non-locality. Additionally, we provide 
a new, simpler, formalism to derive the heat radiation 
(see Eq. ( fl9| ) below). To this end, we start from quan- 
tum thermal fluctuations inside the object following the 
theory of fluctuational electrodynamics introduced over 
60 years ago by Rytov [36 . Then, we derive the heat 
radiation of the experimentally important case of a cylin- 
drical object, discussing polarization effects for different 
conducting and insulating materials, as well as asymp- 
totic limits. This detailed dicussion hence provides much 
deeper understanding of material and size effects com- 
pared to the brief introduction of dielectric cylinders in 
Ref. [28] . We derive the so far unknown heat radiation of 
a cylinder made of uniaxial material, where the symmetry 
axis of the material and the cylinder coincide. We apply 
this formula to introduce a simple model for the heat 
radiation of multi- walled carbon nanotubes (MWCNT). 

The paper is composed as follows: In section|TTJ we give 
the general formalism for heat radiation of arbitrary ob- 
jects and apply it to the cases of a cylindrical object. Ad- 
ditionally, we find the radiation of an anisotropic plate, 
generalizing the known results |36j for the isotropic case. 
Sec. |III| finally gives the specific form of the scattering 
operator for a cylinder made of uniaxial material, needed 
to evaluate the formula of Sec. |II B| In Sec. |IV| we dis- 
cuss the total radiated energy by cylinders made of differ- 
ent materials, such as dielectrics, metals, and MWCNT, 
putting special focus on the degree of polarization of the 
radiation. Finally, in Sec.|Vj we discuss the spectral den- 
sity of the energy emitted by these materials, a quantity 
which might be most accessible in experiments. We close 
with a summary and discussion of our findings in Sec.|VI| 



II. HEAT RADIATION IN TERMS OF 
SCATTERING OPERATOR 



A. General formalism for arbitrary objects 

Consider an object with homogeneous temperature 
T fcj placed in vacuum, enclosed by an environment (e.g. 
a cavity much larger than all other scales in the system) 
at temperature T env . In global equilibrium, i.e. with 
T bj = T env = T, the autocorrelation function C of the 
electric field is related to the imaginary part of the dyadic 
Green's function Gij of the object by the fluctuation- 
dissipation theorem (FDT) [551137]. 



Ct 3 q (T) ee {Ei{u; v)E* {w; v')) eq ee (E(w; r) ® E* ( W ; r')> 



fij 



[a T (u) +a (u})} — ImG tJ (w;r, r'), 



where ® denotes a dyadic product. ax(uj) = 

- y (exp[faj I 'UbT] — 1) _1 describes the thermal con- 
tribution to quantum fluctuations, compare Eq. |T]). 
The zero point fluctuations, which contribute ao(w) = 
w ft(47r) ^ are m( Jependent of the object's temperature and 
do not contribute to heat radiation. Hereafter we use the 
operator notation G = Gjj(u;; r, r'), where operator mul- 
tiplication implies an integration over space as well as a 
3x3 spatial matrix multiplication, e.g. for the operators 
A and B (using Einstein summation convention), 



fe (r,r") = J dV4 y (r,r')fi 3i (r',r"). 



The Green's function is the solution of 

,2 



HO 



V- 



(4) 



(5) 



which follows because the electric field obeys the 
Helmholtz equation, Eq. ^ below. Here, I = VxVx 

describes free space, and V = ^(1 — £+V x ^4 — 1^ Vx) 

is the potential introduced by the object, c and qj are the 
complex (possibly nonlocal) dielectric permittivity and 
magnetic permeability tensors of the object. For isotropic 
and local materials, they reduce to scalars (e.g. c = el). 
Go is the Green's function of free space. Using the iden- 
tities ImG = — GImG _1 G* and ImV = In^G" 1 - G^ 1 ) 
[37], which can be found from Eq. ([5]), we obtain 

c 2 

C eq (T) = C + C(T) - a T H— GImGo 1 G*, (6) 
c 2 

C(T obj ) = -a To6j (w)^GImVG* 



r(w) J d\'d?r"G l3 {v 1 v') 

obj 



xIm^(r',r")G^(r",r"'), 



(7) 



where Co = ao(w)^ImG is the zero point term. Equa- 
tion ^ shows two different finite temperature contribu- 
tions to the electric field in equilibrium. C(T) contains 
an explicit integral over the sources within the object, 
as ImV is only nonzero inside the object, and we iden- 
tify it with the desired heat radiation from the object. 
The expression in Eq. |7| can be shown to be identical 
to expressions in the literature for both complex electric 
and magnetic permeabilities |36[ 137] , where, in general, 
one has two terms, including Ime and Imp, respectively. 
The introduction of the potential V appears useful here, 
as it allows for a compact notation including both terms. 
The third term in Eq. 

r2 



C env (T env ) = -a T _(u J )- 



— GImG 



(8) 



(3) 



is the contribution sourced by the environment. As a 
specific model for environment, consider the objects en- 
closed in a very large black cavity maintained at tem- 
perature T env . This latter identification can be corrob- 
orated on a different route by introducing a cold object 
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into the thermal background field Eo sourced by the en- 
vironment, with field correlator given by (Eo ® Eq) = 

fl T„„ (w)fjImGo. At this point, it is useful to introduce 
the T operator or scattering amplitude T[35J HU] of the 
object. It relates the homogeneous solution (also some- 
times called the exciting field [40]) E h (for V = 0) of the 
Hclmholtz equation 



+ V- 



E = 0, 



(9) 



to its (inhomogeneous) solution Ej/, with the object 
present. This solution can be stated in terms of the 
Lippmann-Schwinger equation 



E,,, = (1 - GoT)E,, 



(10) 



With this equation, the above introduction of the cold 
object into the free environment field is readily done and 
C env is then the correlator of the field E s , the inhomo- 
geneous solution with the object present, 

C env (T env ) = (E s ® EJ) = (1 - GoT)(E ® E* ) 

c 2 

x (-TG* + 1) = -a Te _(w)^GImG ( 7 1 G*, (11) 

in agreement with Eq. ^ . Here we used the identity [35J 
G = G - G TG . (12) 



Equation (11) highlights the physical interpretation of 
(jenv . j-j- j g ra( Ji a tion sourced by the environment and 
scattered by the object. 

Having found the contributions of the different sources 
(environment and object), one can now vary the tem- 
perature of these independently in order to arrive at the 
field outside the object when its temperature is different 
from that of the environment. If T env = this field cor- 
responds to the heat radiation of the object. To this end, 
we notice that it is not necessary to derive all the terms 
in Eq. ^ as explained in the following: The explicit 
expression for C(T fy) in Eq. ^ contains the Green's 
function with one argument inside and one argument out- 
side the object. While this function can be in principle 
derived, we find it more convenient to express C(T t,j) 
in terms of the Green's function with both arguments 
outside the object, as it is directly linked to the scat- 
tering operator by Eq. (12 1. Therefore, it is interesting 
to note that C env has all the sources outside the object 
and hence can be found in terms of this Green's function. 



While this is already obvious in Eq. ( 11 ) we additionally 
present a more rigorous way to derive C env . The envi- 
ronment sources, described by e env , can be thought of as 
being everywhere in the infinite space complementary to 
the object, infinitesimal in strength (environment "dust" 
[37]), i.e. e env -» 1. C env in Eq. ^ can hence be written 

C e ™(T em) )=a Te „» 

lim / d 3 r'G ik (r,r')JmE env G* jk (v",r'), (13) 

£e ™~ 5-1 ./outside 



which is identical to Eq. (4) in Ref. [28^-Here, we intro- 
duced a Green's function G with V inside the object and 
£ e nv outside. This is a simple modification of G as a fi- 
nite s env — 1 only changes the external speed of light so 
that c in G is replaced by c/yje env . 

Finally the heat radiation of the object at temperature 
T a bj can now be found by solving Eq. ^ for C(T t,j), 



C{T obj ) = fflro6 »-ImG - C env (T obj ), 



(14) 



where G is found using Eq. JT2j . Note that C env (T) can 
be derived from either Eq. ( 13 1 or directly from Eq. (|11|). 



For the case of the cylinder, we present below the former 
derivation in detail, and briefly sketch the latter starting 
from Eq. (fTTj) . 



We emphasize again that the field emitted by the ob- 
ject in Eq. (14) is fully expressed in terms of the Green's 



function with both arguments outside the object. In case 
one is only interested in the total heat emitted, the first 
term, i.e., the equilibrium field need not be derived, as 
it contains no Poynting vector. If, on the other hand, 
the radiation of the object is scattered at other objects, 
e.g. in order to compute heat transfer or nonequilibrium 
Casimir interactions, the full expression ( 14 ) has to be 
kept. 



B. Heat radiation of a cylindrical object 



In order to compute the heat radiation of a cylindrical 



object (denoted by subscripts c), we apply Eq. ( 14 ), eval- 
uating the environment contribution by use of Eq. ( |13[ ). 
Afterwards in this subsection, we briefly sketch the 
derivation via Eq. ( |11[ ) . In the cylindrical geometry (with 
parallel, radial and angular coordinates z, r and 4>, re- 
spectively), the free Green's function Go is expanded 
in cylindrical vector waves, RM„ ^ and RN,, ^ corre- 
sponding to Af-polarized and TV-polarized regular waves 
[4*0] . see App. A. These are indexed by the multipole or- 
der n and k z , the component of the wavenumber k = lj/c 
along the cylindrical axis. For outgoing waves we use 
~bJL n ,k z and N n ,kz accordingly. In this basis, the T oper- 
ator of a cylindrical object is diagonal in n and k z , but 
couples different polarizations. Its entry T^ k p relates 
the amplitude of a scattered wave of polarization P' in 
response to an incoming wave of unit amplitude and po- 
larization P, with P, P' £ {M, N}. More precisely, the 
application of the T operator in Eq. ( 10 1 on regular cylin- 
drical functions reads, 



JoT c RP„^=^T^P; 



(15) 



With these definitions and Eq. (A.2), the Green's func- 



tion of the cylinder is easily found, by use of Eq. (12), 
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as 



Jo - G T C G = Go + ]T E (-1) 



p ' pi n——oo 



t\rnPP' 



(16) 



8,7 



-P„ >fc >)®P' 



When performing the integration in Eq. (13), we note 



that Go(r,r') in Eq. (A. 2) is separated into two pieces, 
corresponding to r < r and r' < r. The contribution 
of a finite region vanishes asymptotically in the limit of 
Senv — > 1 and can thus be neglected without changing 
the result. We can hence restrict the integration range 
to r' > r, r", where we have to use exclusively one of 
the pieces. In general, one can restrict the integration in 



Eq. (13) to £(r') > £(r),£(r"), where £ is the component 
which distinguishes the two expansions of Go- 

Due to the orthogonality of two basis sets of the wave 
functions, the integrations over polar angle <fi' and cylin- 



drical axis z' yield 27rJ„ i „, and 2irS(k z — k' z ), respectively, 
and we are left with only one term for each polarization 



in Eq. (13), 



lim / r'dr' 



P-n,-fc. (r') 



2 2c 2 1 



(17) 



where P nj t, has analogous form to P n .k z with the 
wavenumber ^/e env uj/c instead of uj/c. Also, ". . . " repre- 
sent higher order terms in e env — 1. This equation holds 
for k\ < w 2 /c 2 , for k 2 > uj 2 /c 2 all terms are of order e° env 
and do not contribute in Eq. (13), a manifestation of 



the fact that the environment radiation does not contain 
evanescent waves. 

The radiation from the environment after scattering at 
the cylinder then reads 



Cr (T en „)(r,r") = a Te „» EE/ ^W*) ® HP^r")^' + ^.W ® K^")<C 

p t p> n=—oo 

+ P„, fci (r) ® RP- K (r")TX + P»A (r) ® P^^ (r") E ^T^f '* ) - 

(18) 

I 



where P stands for the polarization opposite to P. Phys- 
ically, Eq. ( [18] ) describes the thermal field for the case of 
an environment at temperature T env and a cold cylinder 
(T c = 0). It can also be derived via Eq. (11), by noting 



that the radiation of the environment without cylinder 
present can be given in closed form, 



E (g> Eq) = a Te „, (w) -J Im G = a Tenv (u) E E 



8-7T 



P n=— oo 

rp: 



n,k z 



(19) 



Application of (1 — GoT c ) from both sides (compare 
Eq. (11)), using Eq. (15) immediately leads to Eq. (18). 



This simple route towards the environment radiation 
(and hence the radiation of the object) has not been pre- 
sented before. 

Since we are only interested in the energy emitted by 
the cylinder, we do not explicitly derive the equilibrium 
field in Eq. (14) as it contains no Poynting vector. Equa- 



tion (14) thus states Kirchhoff's law, that the energy ab- 

T, T c = is the 
: 0, T c = T. 



sorbed by the cylinder in the case T el 
same as the energy radiated by it for T ei 
This is a special case of detailed balance in equilibrium, 



which generally states that the absorption coefficient of 
an object equals its emission coefficient. The formalism 
described here hence also provides a convenient route to 
find the absorption coefficient of arbitrary objects. Tech- 
nically, due to these considerations, the Poynting vector 



(S(r)) 



c 

4tt 



Re[(E(w,r)xH*(w,r))] 1 (20) 



of the field in Eq. ( 18 ) gives complete information about 



the net energy flux for any temperature combinations. It 
can be derived via B(w, r) = x E(w,r) as well as 

relation (A. 3). The power \H C \ radiated per length L of 
the infinite cylinder in the general case of finite T env and 
T c is finally given by [2"8] . 



\Hc 
L 



uidu) 



1 



1 



E E 

P=M,N n=-oo 



(21) 



Obviously, if T c < T env , there is a net energy flux into 
the cylinder. In the following, we consider exclusively the 
case T env = (and denote T c = T), for which case the 
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energy flux is referred to as heat radiation of the cylin- 
der. The expression ( 21 ) can be split into two terms each 



representing a different polarization of the corresponding 
electric field. Specifically, the term which describes polar- 
ization parallel to the cylinder is given by the ./V-modes, 



\H N \ = |g||] 
L L 

pUJ I C 



E 



n— — oo 



—U) I c 



oo 

" e 

NN-i , 



ujduj 

k B T _ I 



l«| 2 



(22) 



\rpNM \2 \ 



whereas the term responsible for the polarization perpen- 
dicular to the cylindrical axis is given by the M-modes, 



M 



oo 

E 

n=— oo 



L 



ujcLo 



e K B j 



dk z {MT^] + \TXV + \ T nZ?) 



(23) 



In the following we use a standard definition of the degree 
of polarization / in order to quantify polarization effects, 



N 



\H 



M\ 



H N \ + Iff 



(24) 



In case one prefers a description in terms of the scatter- 
ing matrix S [Ml EE] with S^T = 2Tf' p + S P )P >, the 



radiation in Eq. (21 1 can equivalently be written 

1 



L 



h 

oo 



ojduj 



1 



g Kg 1 



- 1 



E 

P^pi n— — oo 



E 



UJ / c 



dk z (\sX\ 2 -s 



p,P' 



(25) 



Furthermore we note that the result for the perpendic- 
ular emission of the cylinder given in Ref. |24) can be 
recovered from Eq. (21) by restricting to k z — 0, taking 



into account waves normal to the cylindrical axis only. 
In this case the T operator is diagonal in polarization P. 



Limit of large radius (radiation of a plate of 
anisotropic material) 



Casimir forces between a uniaxial plate and single-walled 
carbon nanotubes [42], but we have not come across an 
explicit result for radiation of a plate. For materials with 
anisotropic electric or magnetic response, the Fresnel co- 
efficients are not diagonal in polarization s and p (see 
Ref. [43] for these lengthy coefficients), but take the gen- 
eral form Q for a scattered wave of polarization Q' in 
response to an incoming wave of polarization Q. Note, 
s(p) polarization corresponds to the wave with the elec- 
tric field vector perpendicular (parallel) to the plane of 
incidence. Thus, the heat radiated per surface area (the 
Poynting vector S) can in this general case easily be 
found from Eq. ( |14[ ), where, using a plane- waves basis 
|40) . the steps are similar to the ones performed for the 



cylinder (Eqs. ( 16 )-( IS )) 



S = 



8tt3 Jo 



duj- 



d 2 k ± 



kj_ <UJ I C 



j _ i \ r Q\2 j_ |„QQ|2 



(26) 



Here, k± is the wave-vector component parallel to the 
plate. For r®Q = 0, this equation reduces to the well- 
known one for isotropic materials |36j . The expression 



( 26 1 can be rewritten in terms of M and N polarization 
for cylindrical geometry. If we define 4> to be the angle 
between the optical axis (which, in order to describe the 
radiation of a thick cylinder in terms of the one of a plate, 
is parallel to the plate surface) and the intersection line 
between the plane of incidence and plate itself, then we 
can write the limiting values of the M and N components 
of the cylinder radiation as, 



Sm 



h 
1 

2 



dui — 
e 

„s|2 



27T /* aJ / c 

dtp / k±dk± 
' 7 - l Jo Jo 



i-^{(k1 2 + k sp l 2 )cos 2 </> + (|7 



|r ps | 2 ) sin 2 . 



(27) 



h f°° uj f 27r f u/c 

S N = —. j / duj— ^ / d(f> / k ± dk± 

Stt Jo e 1 ^ - 1 Jo Jo 



1 - o {(l r '| 2 + l r " P | 2 ) sin2 + rP ! 2 + \ rPS \ 2 ) cos2 ' 



(28) 



For large radius, the radiation of a cylinder is asymp- 
totically identical to the radiation of a plate (a semi- 
infinite planar object) of same surface area [35]. For the 
case of a plate made of isotropic material, the heat ra- 
diation is well-known |36| . Nevertheless, we will below 
study the heat radiation of a cylinder made of uniaxial 
material (as a simple model for carbon nanotubes), see 
Fig. [T] in which case the limit of large radius is a plate of 
uniaxial material. Recent literature discusses heat trans- 
fer between plates of uniaxial materials [JT] as well as 



Note r sp = r ps when the optical axis is parallel to the 
plate. 



III. 



T OPERATOR FOR A CYLINDER MADE 
OF UNIAXIAL MATERIAL 



In section [IIB[ we derived the heat radiation of a cylin- 
drical object expressed in terms of its T operator which 
is known analytically |1U FH] . One aim of this paper is 
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FIG. 1: (color online) Infinitely long cylinder made of a uni- 
axial material. The symmetry axis of the cylinder coincides 
with the symmetry axis of the diagonal dielectric tensor (the 
z-axis). Electromagnetic waves radiated by the cylinder are 
denoted by M and N for Af-polarized (perpendicular) and N- 
polarized (parallel) respectively. Note, generally M-polarized 
waves have components along both azimuthal and radial di- 
rections. Also, despite iV-polarized waves have non-zero com- 
ponents along all three basis directions, it is only the com- 
ponent parallel to the cylindrical axis to contribute to the 
Poynting vector. 



to study the radiation of a cylinder of uniaxial materia l 

IVC] ). 



as a simple model for carbon nanotubes (see Sec. 
The corresponding T operator, which seems unavailable 
in the literature, will be derived here for the case depicted 
in Figjl] where the cylindrical axis coincides with the op- 
tical axis. This is done by solving the scattering problem 
in Eq. ( jlOj ), which amounts to satisfying the boundary 
conditions for the electromagnetic field at the cylinder's 
surface. The setup is described by isotropic local mag- 
netic permeability fi(u>) and the following tensor for the 
anisotropic, but homogeneous and local dielectric func- 
tion inside the cylinder, 



e r (w) 
e{uj) = ( e r (u) 
e,(w) 



(29) 



The electric displacement field inside the cylinder can 
then be written as 

D = e r (u) (E r e r + E^) + e z (oj)E z e z , 

where e r , and e z correspond to unit vectors in cylin- 
drical coordinates. 

Importantly, uniaxial materials split the incident beam 
into ordinary and extraordinary ones [45] . In our geom- 
etry, cylindrical M-polarized waves correspond to ordi- 
nary ones and propagate according to the dielectric func- 
tion e r . The ./V modes correspond to extraordinary waves 
and propagate according to an effective dielectric func- 
tion which depends on the direction of propagation. 

The resulting T operator components T^jf , as defined 
in Eq. (15), are expressed in terms of Bessel functions, 



J„, and Hankel functions of first kind, (see App. |B 
for the detailed derivation), 



J n (qR) A1A4 
H^\qR) AiA 2 



K 2 



(30) 



rpNN 

1 n,k z 



rpNM _ rpMN _ 

n,k z n,k z 



J n {qR) A 2 A 3 -K 2 



(31) 



where 



2i K 1 

TT^(qR) 2 [H^iqR)} 2 AiA 2 - K 2 ' 

(32) 

(33) 



J'nilNR) 



q N RJ n {q N R) £zqRH^\qR) 



1 Hi lY (qR) 



Hi lY (qR) 



q M RJn(qMR) ^qRH ( n\qR) 



A, 



J'nilNR) 



1 J' n {qR) 



q N RJ n (q N R) E z qRJ n (qR) 



J'tMmR) 



1 J' n (qR) 



q M RJ n {qMR) fJ,qRJ n (qR)' 



and 



K = 



(34) 



(35) 



(36) 



(37) 



k = uj/c and q — \Jk 2 — k 2 are the wave- vector mag- 
nitude in vacuum and its component perpendicular to 
the z-axis respectively. qM = y/ e r fj,k 2 — k 2 and q^ = 
\J £ z j E r \J E r \ik 2 — k 2 are the wave-vector components 
perpendicular to the z-axis for the M-polarized ordinary 
and A^-polarized extraordinary waves inside the cylinder 
respectively. 

As required by continuity, the above T matrix can 
be easily reduced to the isotropic case when e r — 
e z = e. Then the expressions simplify to qM = 
qN = \J s\ik 2 — k 2 = q 1 and K = nk z c(l/q' 2 — 

l/9 2 )/(v^^ 2w )' an< ^ our resurts reduce to the known 
forms for an isotropic cylinder [TU [II] . 



IV. EXAMPLES AND ASYMPTOTIC RESULTS 

In this section we numerically, and analytically, study 
the radiation of a cylinder for different material classes. 
We start with isotropic dielectrics and conductors and 
finally present the case of a uniaxial material using the 
"in-layer" and "inter-layer" response of graphite to model 
MWCNT. We consider fi = 1 for all studied materials. 



A. Dielectric cylinder 

Figure [2] illustrates the result of the numerical calcu- 
lation of the radiation by Si02 and SiC cylinders for 
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FIG. 2: (color online) The normalized heat radiation as a 
function of radius R for (a) Si02 and (b) SiC cylinders at 
temperature T = 300K. Calculations were performed using 
Eq. ( |21[ |, and analytical expansions, Eq. (38 1 and Eq. (39 1. 
The horizontal lines show the radiation of (a) Si02 and (b) 
SiC plates. At and the smallest skin depths 5 in the relevant 
frequency range are marked on i?-axis. 



T = 30QK normalized by the Stcfan-Boltzmann value, 
Eq. (2|). 

For SiC>2 optical data was used, whereas for SiC the 
following dielectric function was taken [46] . 



where e M = 6.7, uj lo = 0.12eV, uj T o = 0.098eV", 7 = 
5.88 x 10- 4 eV. 

We see the effects discussed in Ref. [28 involving the 
three length scales: radius R, thermal wavelength At, 
and skin depth S — c/ (Jm^/euj), where the latter de- 
pends on frequency. If R is the smallest scale, i.e., much 
smaller than the smallest relevant skin depth and At, the 
radiation is proportional to the volume of the cylinder. 
In this case, radiation emitted inside the cylinder will 
hardly be reabsorbed on its way out so that all regions 
of the cylinder contribute equally to the emission. The 



other asymptotic behavior is approached when R is the 
largest scale, i.e. much larger than the largest relevant 
skin depth and At- Then only the cylinder surface con- 
tributes to the radiation which approaches the values of 
a plate of equal surface area as seen in the figure. 

The length scale At sets (via the Boltzmann factor) 
the range of relevant wavelengths of emission. Neverthe- 
less, for dielectrics there is a fine structure to this range 
given by the resonances of the material. In case of R <C 5, 
the cylinder emits predominantly at the resonance wave- 
lengths of the material, where 5 is minimal. On the other 
hand, for R ^> {S, At}, the emissivity is strongest in re- 
gions where e f=a 1 (compare the plate emissivity). 

In general, one might expect resonance effects when 
the emitted wavelength is of the order of Rj (2-7r) (similar 
to Mie resonances for a sphere [TT1 IT2"]). Due to the 
contribution of all wavelengths, these are smeared out in 
the total heat emitted. For SiC>2 in Fig. [2] the fact that 
the emissivity exceeds the plate result for R w 20/im 
might be connected to such resonances. 

Figure [2] shows imprints of the dielectric function of 
SiC which has a sharp strong resonance (leading to a 
small 5), but apart from the resonance, SiC is almost 
black (i.e. has very large S) in our frequency range. The 
two regimes R <C {5, At} and R 3> {S, At}, where the 
cylinder radiation follows the discussed asymptotic laws 
are far apart, and very large radii R are necessary to ap- 
proach the classical plate result. In Fig. [2j the radiation 
for R = 150/im is still distinctly different from the asymp- 
totic result. This might be advantageous for experiments, 
as a SiC cylinder does not have to be very thin in order 
to observe deviations from the Stefan-Boltzmann law. 

In the limit of R <C {5, At}, the radiation can be stud- 
ied analytically, see App.[Cj where we find the following 
asymptotic laws for the two polarizations, 



ton M 

-R«{>5,A T } A 



_ 1 

~ 6 
x Im 



W) 

e(u>) 2 + 2e(w) 
e{u) + 1 



A a T (u)cR 



(38) 



ton ™ 

i?<S{5,A T } A 



dw 



ar(w)ci?Im 



e(ui) + 1 ' 



(39) 



Note the linear increase with R, corresponding to the pro- 
portionality of the unnormalized radiation to the volume. 
The numerical evaluation of these asymptotic forms has 
been added to the graphs in Fig. [2j where the agreement 
for small R to the full results is visible. As expected 
and seen from the equations above, cylinders with purely 
real dielectric functions (or more general with real po- 
tential V) will not radiate, because the dissipative prop- 
erties of the material are responsible for the heat radia- 
tion (in accord with FDT). This holds for any R. The 
iV-polarization given in Eq. ( 38 1 dominates over the M- 
polarization in Eq. ( p9| if 

(40) 



[Ree(w) + l] 2 + [Ime(w)] 2 > Ime(u;) 



holds, which is the case for most materials. Further in- 
sight can be gained by additionally requiring the tem- 
perature to be so low that one can expand the dielectric 
function, i.e., At 3> Ao, where Ao is the wavelength of 
the lowest resonance of the material. In this case 1551. 



s{u) = £ + 



I h 0(L) ), 



(41) 



with e an d Kn real. Plugging (41) into Eqs.(38) and 



(391 the frequency integration can be done and we have, 



lim 

i?<S{5,A T },A <A T 



\H 



4tt 4 hc 2 X in R 
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(eo + If 



lim 

fl«{<5,A T },Ao<A T 



M 



8?r 4 Hc 2 X in R 1 



A 



(»:-! 



A/J1 



(eo + 1) 5 



(42) 



(43) 



Interestingly, to lowest order in T, both components scale 
as T 6 and hence fundamentally different from the Stefan- 
Boltzmann law in Eq. pi), which scales as T 4 . 

The degree of polarization / in Eq. (24 1 is finally given 

by 



lim / 

fl<{t?,A T },A <A T 



4 



2e -3 



si + 2e Q + 9 



(44) 



Although the condition A <C At is not fulfilled in Fig. [2] 
Eq. (44) still gives a hint to why SiC has higher degree 



of polarization compared to Si02 , as £n is tangibly larger 
for SiC and Eq. ( 44 ) monotonically increases with £n (for 
so > 0). 



B. Well conducting cylinder 

Conductors differ from insulators by a significantly 
smaller skin depth 5, leading to very different radiation 
characteristics. In this section, we first study the radia- 
tion using a simple Drude model with the parameters of 
gold in order to highlight the different limiting behaviors. 
Then, we turn to tungsten which has been extensively 
used in experiments [16] [24] . 

The Drude model for gold [47], 



£Au(w) 



uj(uj + iuj T ) 



(45) 



1, LU V 



9.03eV and w T 



has the parameters 
2.67 x 10- 2 eV. 

Figure [3] shows the numerical result for the total radi- 
ation by a gold cylinder at 300K and 30iC We observe 
a behavior drastically different from the ones in Fig. [2] 
the radiation, normalized by surface area, can be many 
orders larger than expected from the Stefan-Boltzmann 
law. More precisely, it increases with decreasing R, has 
a maximum at R w S (for Drude model of gold ( 45 1 , 
6(oj) has no sharp minimum, and we show the skin depth 
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FIG. 3: (color online) The normalized heat radiation as a 
function of radius R for Au cylinder at temperatures (a) 
T — 300A" and (b) T = 30K. Calculations were performed 
using the Drude model [47] with Au optical properties. Full 
numerics, Eq. \21\ , analytical expansion, Eq. ( |46| l, and the 
approximation of Eq. (47 1 were used. Black horizontal lines 



indicate the radiation of Au plates at corresponding temper- 
atures. At and the skin depth 5 in the relevant frequency 
range are marked on the _R-axis. In the insets the degree of 
polarizations, |24|, are plotted. 



corresponding roughly to thermal wavelength) and ap- 
proaches the laws ( 38 ) and ( 39 ) only for very small (un- 
physical) radii. For R ^> \ T , the result of a gold plate 
is approached. The large values of radiated power, com- 
pared to the Stefan-Boltzmann law, can be explained by 
the interplay of two effects: a large imaginary part of 
the dielectric function gives rise to strong radiation from 
every volume element of the conductor. On the other 
hand, it also leads to a very small skin depth such that 
most radiation is reabsorbed inside the conductor, and 
only a thin surface layer contributes to the radiation. In 
the region where R rj 5, one has very strong emitting 
elements which all contribute to the total radiation, and 
hence emission normalized by surface area is maximal. 
Interestingly, when the conductivity goes to infinity, the 
effect of vanishing skin depth is stronger than the effect 
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of increasing radiation such that the emissivity vanishes 
as l/\/e. In this case, the classical (plate) limit goes to 
zero and the maximum in the curve shifts to smaller and 
smaller R. 



The insets show the degree of polarization (24 1 as a 
< \m thp radiation from the 
For R > At, the 



function of R. For R < At, the 
gold cylinder is fully ./V-polarized. ror n ^ at 
degree of polarization becomes negative and asymptot- 
ically approaches zero for R/Xt — > oo. These qualitative 
features agree with experimental [TrjHTOl |2"4"] as well as 
theoretical! 24 studies on radiation of thin wires. 



FIG. 4: (color online) The normalized heat radiation as a 
function of radius R for tungsten cylinder at temperatures (a) 
T = 298K and (b) T = 2400K. Calculations were performed 
using Eq.(21| with dielectric function (48 1. Black horizontal 



lines indicate the radiation of tungsten plates at correspond- 
ing temperatures. At and the skin depth 8 in the relevant 
frequency range are marked on i?-axis. In the i nset s the de- 
gree of polarizations are given using expression (24 1. 



For conductors, the appropriate limit for an analytic 
expansion is At ^> R 3> S (where S is of the order of 
nanometers for good conductors). In this case, the lead- 
ing element of the T operator is T^^, see App. pi The 
resulting radiation is then completely iV-polarizea and 
reads 



lim 



\H r 



or(w) 



(2tt) 2 u; 



tt/2 



d6- 



2Re[l/ve] cos 3 9 



|cos 2 9(2^ E - iit)ujR/c - 2ij^e + 2 cos 2 9uR/c log [cos 9ujR/2c 



(46) 



where je ~ 0.577 is Euler-Mascheroni constant and 9 is 
the angle of incidence so that k z — ksva.9. 

As described above, we see here explicitly that the ra- 
diation (in the considered range of radii) vanishes for 
|e| — > oo. This expression cannot be further simplified as 
the integrand diverges at 9 = tt/2 if we omit the small 
term — 2i/y/e in the denominator. A result almost simi- 
lar to ( 46 1 was obtained by Rytov [48] (whose derivation 
was restricted to Ree = and line 3> 0), which differs 
by the absence of the first term in denominator. We em- 



phasize however that this term is necessary for accurate 
results in the considered limit. 



Rytov suggests a further simplification of his expres- 
sion by setting cosf5 = 1 in log [cos 9tdR/2c], so that 
integration can be performed analytically. Omitting 
cos 2 9(2~/E—iTr)LoR/c in denominator, the following result 



10 



for radiation can be obtained after integration, 

\ha r 00 du J' 2 



lim 



A r ».R»<5 A 



o (2tt) 2 2(2tt) 2 w 5 / 2 



x a T (w)Re[l/v / e]- 



11/4 



i?log[a;i?/2c]| 3 /2- 



(47) 



Note that Eq. ( |47[ ) was derived for any well-conducting 
(non magnetic) media, whereas the corresponding result, 
Eq. (IV. 35) in Ref. [J5], only considers the case of a 
purely imaginary dielectric function (and fi =/= 1). Due to 
this the integrand of Eq. (47 1 differs from Eq. (IV. 35) 
in Ref. 05] by Re[l/^}\e\^ 4 instead of (/Valine) 1 / 4 
(which agree for Ree = and fj, = 1). Also, we empha- 
size that Eq. |47] ) is an approximation whereas Eq. ( 46 ) 
is exact asymptotic limit for At 3> i? S> S. 

Figure [3] provides a test of these approximations, 
demonstrating that for T = 30K, where the ratio of 
V : 



and 8 is larger than at T = 300K, Eq. (|46j) de- 
scribes the full solution over roughly one decade in R. 
For T = 300-ftT, the agreement is not as good as the ratio 
of At and S is too small. Equation (47) gives a rough 
estimate of the overall dependence on R in Fig. [3j but 
the values differ by roughly a factor of 10 from the ex- 
act results. Moreover, above some threshold value of R 
we cannot obtain a finite radiation rate from Eq. (47) 



because of the divergent term due to log of unity in the 
denominator of the integrand. 

After this (more theoretical) discussion of gold, where 
we used the same e for both temperatures in order to 
demonstrate the pure effect of temperature via the Boltz- 
mann factor, we turn to the experimentally relevant ma- 
terial tungsten, at relevant temperatures of T = 298K 
and T — 2400-ftf as shown in Fig. |4j We use the corre- 
sponding dielectric function for tungsten|49j. 



A 2 



Kq p X 



i5pX S pX 



A 2 
2ttc€q 



E 



1 ^rq 



(48) 

where A is the wavelength in vacuum, c is the velocity 
of light and en is the permittivity of vacuum in SI units. 
The remaining parameters are listed in Table [Tj 

The overall radiation of tungsten at 298K is very sim- 
ilar to gold at 300K in Fig. [3j At high temperature 
T = 2400/-T , the increase of the normalized radiation 
over the Stefan-Boltzmann law is reduced. We attribute 
this to a smaller conductivity at T — 2400X. We also 
observe in the insets that the zero in the polarization 
curves shifts by roughly a factor of 10 when comparing 
the two temperatures. This manifests that the zero in 
the polarization curves is mostly a function of R/Xt as 
At is also shifted by almost a factor of 10. Furthermore, 
although the polarization in the inset is indistinguishable 
from unity at small R, we note that the ratio of V and 
M polarizations is a factor of 10 3 larger for T = 298K 
compared to T = 2400 K. 
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TABLE I: Optical data for tungsten from |49| . Temperature 
T is in Kelvins. Conductivities (<ti, etc.) are in units of 10 6 
ohm - m -1 . Wavelengths (A r i> etc.) are in /xm. Kq\, etc. 
and Si, etc. are dimensionless. 
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FIG. 5: (color online) The heat radiation of a MWCNT as 
a function of radius R, normalized by the Stefan-Boltzmann 
result, at T — 300iV Contributions from two polarizations 
are indicated. Corresponding thin curves without boxes rep- 
resent the heat radiation for "isotropic graphite" with dielec- 
tric function (e z + e r )/2. Horizontal lines of different colors 
indicate the graphite plate classical result for correspondingly 
colored curves. At and the smallest skin depth 8 in the rel- 
evant frequency range are marked on the i?-axis. Note that 
the smallest skin depths corresponding to both e r and e z are 
approximately equal and labeled here by 8. In the inset the 



degree of polarization is given using expression ( 24 1 



C. Multi- walled carbon nanotube 

We finally turn to heat radiation of a cylinder made of 
uniaxial material. This can be considered a simple model 
for a MWCNT 50J which is of high importance in modern 
science and technology. As a MWCNT is a wrapped up 
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TABLE II: Optical parameters for the in-layer dielectric func- 
tion of graphite (e z ) from 51 . fj and aj are dimensionless, 
whereas uj t j and Fj are in eV. 
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TABLE III: Optical parameters for the inter-layer dielectric 
function of graphite (e r ) from |51| . fj and a.j are dimension- 
less, whereas cotj and Fj are in eV. 



graphite sheet, we can in a crude approximation regard 
it as a (solid) cylinder described by two different dielec- 
tric response functions: the response along the cylindri- 
cal axis is given by the in-layer properties of the graphite 
sheets, whereas the response perpendicular to this axis 
is approximately given by the inter-layer response. We 
also note that most mineral crystals have uniaxial optical 
properties, and heat emission by these materials might 
open new possibilities for applications |41j. 

Figure [5] shows the heat radiation of a MWCNT for 
T = 300K using expressions (21 1 and (30)-(37), normal- 

oT A A. 



ized as before by the Stefan-Boltzmann law H 
We used the following form for the graphite dielectric 
function|51). 
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n 2 



uj (u) + iT ) 
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3^p 



(49) 



where tt p — VTo^p an d a specific form of 1^ = 

21 

is used, which best describes 



r, exp 



huj—hujj 



the experimental data. For the dielectric function com- 
ponent e z along the cylindrical axis, i.e., the in-layer re- 
sponse, the parameters are uj p — 19eV, T = 0.091eF, 
/o = 0.016. The remaining parameters are given in Table 
[TT) For the dielectric function component e r perpendicu- 
lar to the cylindrical axis, i.e., the inter-layer response, 
the parameters are uj p = 27eV, T = 6.365el^, / = 0.014 
and the remaining parameters are given in Table |III| 

These parameterizations apply to the frequency range 
0.12 — 40eV and 2 — 40e^ for e r and e z respectively, 
but we nevertheless use them for the range of roughly 
0.004 — 0.2eFdue to the lack of simple analytic forms 
for the broader range. Note that At corresponds to a 
frequency of 0.163eV for T = 300K. 

The overall radiation curve is in between the character- 
istic shapes of dielectrics and conductors (compare Figs [2] 



and [3[): t he regime proportional to volume as in Eqs. ( |38| ) 
and (|39[) is visible for small R in contrast to Fig. [3] On 
the other hand, the strong increase over the plate-result 
and over the Stefan-Boltzmann value, characteristic for 
conductors, is visible as well. These features follow from 



the dielectric functions in Eq.(49 ) carrying a smaller con- 
ductivity compared to gold. 

In order to highlight the effect of material anisotropy, 
in Fig. [5] we also show the thin curves without boxes for 
which we use an isotropic dielectric function, given by 
(e z + e r )/2. We see that the influence of anisotropy on 
the radiation is almost negligible at small R, whereas 
at intermediate R it strongly increases the degree of po- 
larization perpendicular to the tube. Interestingly, the 
asymptotic value of / in the inset is different from 
and takes the value —0.297, an effect purely due to the 
anisotropy which is computed using the result for a plate 
of anisotropic material, Eqs. (27) and (28 1. We note that 



a very thick MWCNT might in fact be best described 
by a plate with optical axis perpendicular to the surface. 
Thus, while Fig. [5] gives the correct radiation for a ma- 
terial with the dielectric properties given in Eq. ( 49 ) , it 
probably only describes MWCNT for small R. Since, 
at small R, the polarization is hardly dependent on the 
anisotropy of the material, we conclude that the polar- 
ization effects for MWCNT [21] are mainly an effect of 
cylindrical geometry rather than anisotropy of the mate- 
rial. 



SPECTRAL EMISSIVITY 



In Sec. IV 



we studied the total heat radiation of a 
cylinder made of different materials. While this is of in- 
terest in connection with efficient heating or cooling, an- 
other quantity which can be more appropriate for direct 
comparison to experiments is the spectral emissivity. In 
this section, we discuss the spectral emissivity for cylin- 
ders made of dielectrics (Si0 2 ), conductors (tungsten) 
and uniaxial materials (graphite) for a fixed radius of 
R = 5/xm. 

The spectral emissivity (density) H u is given by the 
integrand of Eq. (21), 
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E E 
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(50) 



o/c 



dk z (Re[T-l] + \T-j:f + \T^ 



We denote by H u m and H^n the correspondingly polar- 
ized components of . 



A. Si02 

Figure [6] illustrates the spectral density for Si02 at 
T = 300K. The unsteady local fine structure reflects 
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FIG. 6: (color online) The spectral density divided by Stefan- 
Boltzmann law as a function of wavelength A for SiCb cylinder 
of radius R = at temperature T — 300K. In the inset 
the spectral degree of polarization is provided. 



the nature of the optical data which has a number of 
smaller peaks. For short wavelengths (high frequen- 
cies) M-polarized radiation mostly dominates, whereas 
for A > 25/im the TV-polarized radiation starts to prevail 
up to the limit of long wavelengths (low frequencies). In 
the inset the spectral degree of polarization as a function 
of wavelength is shown, where 



■N 



M\ 



uN\ 



(51) 



The two large valleys in this curve are due to the two 
dominant resonances of SiC>2. The fact that the reso- 
nances lead to negative valleys in the polarization is a 
feature specific for the radius chosen. For large wave- 
lengths, the spectral degree of polarization approaches a 
constant value, which can be computed easily using ex- 
pressions (38) and (39), 



lim I u = 

{<5,.R}«A 



|e(w)| 2 +2Ree(cj) -3 
|e(w)| 2 +2Ree(w) +9' 



(52) 



Furthermore, if additionally A 3> Aq holds, the dielectric 



function is described by Eq.(41) and the spectral degree 



of polarization is given then by Eq.(44|. 

Note that the spectral degree of polarization is inde- 
pendent of temperature if the dielectric function is inde- 
pendent of temperature. 



B. Tungsten 

Figure [7] shows the spectral density of radiation for 
tungsten cylinders at T = 298K and T = 2400K. The 
shape of the curves is very similar to Planck's classical 
law due to the Bose-Einstein statistics factor in ar(w). 
The curves peak at values slightly larger than the cor- 
responding thermal wavelengths. For short wavelengths 
A/-polarized radiation is stronger than TV-polarized one, 
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FIG. 7: (color online) The spectral density divided by Stefan- 
Boltzmann law as a function of wavelength A for tungsten 
cylinder of radius R — 5/j.m at temperatures (a) T — 298K 
and (b) T — 2400A". In the insets the spectral degrees of 
polarization are provided. 



whereas for A > 2R, TV-polarized radiation dominates, 
strongly suppressing the M-polarized radiation in the 
limit of long wavelengths. This transition of polariza- 
tion is also manifested in the insets where spectral de- 
grees of polarizations are plotted. In the limit of long 
wavelengths these approach unity, as can be justified by 
results of Scc jlV B} where it is shown that the TV com- 
ponent is dominant in the limit Ay > i? > <). On the 
other hand, in the limit of short wavelengths the spectral 
degree of polarization approaches zero as the cylinder ra- 
diates as an isotropic plate. We emphasize again that the 
spectral degree of polarization is independent of temper- 
ature (if e(lj) is). 



C. Graphite 

Figure [8] shows corresponding results for a graphite 
cylinder at T = 3007^, displaying similar behavior as the 
case of conductors in Figj7| The transition point, where 
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FIG. 8: (color online) The spectral density divided by Stefan- 
Boltzmann law as a function of wavelength A for graphite 
cylinder of radius R — 5/xm at temperature T = 300K. In 
the inset the spectral degree of polarization is provided. 



the polarization changes sign is at A w 25^m. Analo- 
gously to conductors, the spectral degree of polarization 
for a graphite cylinder tends to unity in the limit of long 
wavelengths, i.e. spectral density has polarization paral- 
lel to the cylinder. The range of wavelengths 200— 300/im 
shows an unexpected wiggle in the curves, the origin of 
which is unclear. 

We note that 1^ goes to zero for A — > 0, although 
a finite value is approached for R/Xt — > oo in Fig. [5] 
This is due to the fact that both e z and e r tend to 1 
for w —> oo, and the material is asymptotically isotropic. 
Nevertheless, at A « 5/zm, the polarization is very strong 
compared to Fig. [7J an effect which we attribute to the 
anisotropy of the material. 

Another model describing the spectral degree of po- 
larization of MWCNT's is presented in Ref. [26], where 
we note partly common structure to our description aris- 
ing from the expansions of Eqs. (C.1)-(C.5|. For large 



X/R oo, both the experimental data as well as the 
theoretical predictions of Ref. [26] give values for I w close 
to unity, in agreement to Fig. [8] 



D. Comparing material classes 

Finally, Fig. [| compares the spectral polarization for 
the different materials discussed, where we used simpli- 
fied dielectric functions for dielectrics and uniaxial mate- 
rials in order to illuminate the pure influence of X/R. 

Figure [9] manifests that the overall dependence of the 
polarization on wavelength is quite universal following 
similar curves for all cases shown. Fundamental differ- 
ences are seen in the limiting cases. In the limit of large 
wavelengths all conductors approach unity, whereas di- 
electrics go to a constant value different from unity when 
A > {S,R} (see Eq. ([52])). 

In the opposite limit of small wavelengths, Fig. [^mani- 
fests strong dependence of the polarization on uniaxiality. 



FIG. 9: (color online) Spectral degree of polarization as a 
function of wavelength A for cylinders made of isotropic di- 
electric, isotropic conducting, and anisotropic conducting ma- 
terials for R — 5/im. Black horizontal and orange lines cor- 
respond to long wavelengths limiting values for conductors 
(isotropic and uniaxial) and the dielectric (using Eq. (52l), 
respectively. 



We emphasize that by showing additionally artificial ma- 
terials with very strong anisotropy (factor 199 between 
e r and s z ). Varying this factor, a range from roughly 
0.4 to -0.9 in polarization at A = 10/im can be sweeped. 
Figure [9] clearly illustrates that the smaller e r /e z , the 
smaller the spectral polarization. 



VI. DISCUSSION 

In this article, we derived the general formalism to 
study the heat radiation of a single object held at uni- 
form temperature. While the presented formalism can 
describe the radiation of arbitrary objects, we focus on 
cylindrical objects, providing explicitly the heat radiation 
formula (Eq. pi] )) in terms of its T operator. In order 
to work out the radiation of multi-walled carbon nan- 
otubes in a simple model, we derive the T operator for a 
cylinder made of uniaxial material. We study the cases 
of dielectrics, conductors and MWCNT numerically, dis- 
cussing certain limits (e.g. small radius) analytically. To 
lowest order in temperature T, the heat emitted by di- 
electric cylinders is proportional to T 6 in contrast to T 4 
in the Stefan-Boltzmann law, Eq. ([2|. 

For all materials, the limiting value of very large ra- 
dius is given by the radiation of a plate of equal material. 
In the limit of small radius, dielectrics emit proportional 
to their volume, whereas this regime is not reached for 
conductors in the physically relevant range of radii. In- 
stead, the radiation normalized per surface area can be 
very large for conductors and we observe values as high 
as almost a hundred times the value expected from the 
Stefan-Boltzmann law. The maximum occurs where the 
radius is roughly equal to the skin depth (e.g. a few nm 
for gold) . We note that the validity of using a continuous 
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dielectric function at the scale of nanometers is question- 
able and has to be addressed in future studies. 

The heat radiated by a long cylinder is polarized. All 
studied materials show common overall features of their 
degree of polarization: cylinders with radius smaller than 
the thermal wavelength emit predominantly radiation po- 
larized parallel to their axis, whereas for radii larger than 
the thermal wavelength, the polarization points into the 
direction perpendicular to the cylinder. The limiting 
value of the degree of polarization for small radii is dif- 
ferent for insulators and conductors: while the former 
approach a value less than unity, thin conductors emit 
completely polarized radiation and the degree of polar- 
ization approaches unity. 

The effect of uniaxiality on the polarization is most 
dominant for large radii, where the degree of polariza- 
tion can be tuned over a wide range by adjusting the 
ratio of parallel and perpendicular responses. For small 
radii, uniaxiality hardly influences the polarization. As 
MWCNT fall into the latter size regime, we conclude that 
the polarization measured experimentally for nanotubes 
are mostly a consequence of geometry, and not so much 
of material anisotropy. 

Detailed comparison to experimental data (see e.g. 
Refs. [25H27] for MWCNTs) will be left for future work. 
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Appendix A: Cylindrical harmonics and Free 
Green's function in cylindrical basis 

According to Ref. [40], the cylindrical harmonics can 
be written as, 



RM„, tl (r) 



in 



— J n (qr)e r - J' n (qr)e 4 
qr 



ik. 



ik z z-\-in<p 



RN„ !fe ,(r) = - ik z J' n (qr)e r - : ^J„(qr)e ( ^ 

w [ qr (A.l) 

+qj n (qr)e z ]e^ z+m t 



the cylindrical z-axis respectively satisfying the relation 
q = \Jk 2 — fcj, k = oj/c. J' n corresponds to the first 
derivative with respect to the argument. Furthermore, 
we denote the corresponding outgoing waves by M„ ^ 
and N ?l fe z , which differ from regular ones by replacing 

J„ with the Hankel function of the first kind H„ . 

The above solutions correspond to transverse waves, 
i.e. V • RM n , kz = V • RN, 4j = 0. Moreover, they obey 
useful relations RM n ^ = x RN n ,fc,, RN ni ^ = 
-V x RM„ : fc 2 . These relations are also valid for outgoing 
waves. 

The free Green's function in cylindrical waves 
reads g0], 



Jo(r,r')= E 



P=M,N n= 



idk z 
~8tT 



[RP„ jfcjs (r)(8)P_„_ fea (r')], r'>r 



(A.2) 



[Pni 2 W®RP- 



«.(*')] 



r < r. 



The following relation for propagating cylindrical waves 
is useful for deriving the Poynting vector, 



Re 



# (RP^(r)xVxP^(r)T n p ;f 

+P„, fea (r)x VxRP;* fe (r)Tj; 
= 4<5p P ,Re [T^l] > 



Re 



■ #(P n , fe2 (r)xVxP^ fe ,(r))-e r 
Jo 



4<5p,p'. 
(A.3) 



Appendix B: Scattering of electromagnetic waves 
from uniaxial cylindrical objects 

Consider a wave propagating in an anisotropic 
medium, with dielectric permittivity tensor (29) and 
magnetic permeability n(u>). 

Considering the system's uniaxial symmetry we look 
for wave-solutions in the form of cylindrical harmonics 



(A.li, 



RM^(r) 



x e 



1M T 

ik z z-\-ind 



-J n (qMr)e r - J' n {q M r)e^ 



ik z J' n {q N r)e r -J n (qNr)e^ 

qNr 



+lNqNJn(q N r)e z ] e ik **+ in *, 



(B.l) 



where J n is the Bessel function of order n. RM nj t z and 
RN„ i t i correspond to regular magnetic multipole (TE) 
and electric multipole (TM) waves respectively. Also, k z 
and q are the wavevectors parallel and perpendicular to 



where 7^ is some constant which modifies the N- 
polarized cylindrical harmonic due to uniaxial anisotropy. 
Importantly, we do not care about keeping the harmonics 
(B.l) normalized as that does not influence calculations 
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of T matrix elements in which we are interested. Note 
that the time dependence exp(— iuit) is omitted here. 
The following Maxwell's equations must be satisfied, 



VxVxBCJr) 



V-[eRM™»] =0, 



1 / d 2 



(B.2) 



with analogous relations for RN^ n fc (r) 



Substituting expressions (B.l ) into equations (B.2), we 



obtain the following unnormalized wave solutions, 



Using the definition of the T matrix, we describe the 
scattering process of a regular magnetic wave by the field 
outside the cylinder, which is 



E 



M,out 



, = RM n , fc , + T„X M n ^ + l££ N n , tj (B.5) 



and the field inside the cylinder, 

xM.in AMMjfK/run , a A'M r> \ 



E 



n,k z 



(B.6) 



Analogously, for an incident electric (TM) multipole 
field, the field outside the cylinder becomes 



-Jn(qAir)e r - J' n (q M r)e^, 



ik z z-\-in(p 



RN™, (r) 



ik z J' n (q N r)e r - J n (gjvr)e^, 



+—qNJ n (qNr)e 2 



(B.3) 



ik z z+in<j) 



where qM and q^ are the wave-vector components per- 
pendicular to the z-axis for the two solutions respectively, 



qM 



\/s r fj,k 2 - k% q N = \Js z /e r \/e r nk 2 - fcf. 

(B.4) 

One solution is an ordinary wave and is M-polarizcd, 
whereas the other one is called an extraordinary wave 
and possesses TV-polarization [45] . 

In order to solve the scattering problem for the cylin- 
der, we expand the electromagnetic field in cylindrical 
basis (A.l) and (B.3), outside and inside the cylinder 



respectively. The expansion coefficients for the field in- 
side and outside can be obtained by matching boundary 
conditions at the cylinder's surface for field components 
tangential to the surface. 



E 



N,out 



RN 



n,k z 



T™ k N M n , k ; 



T NN ~N , 



and the field inside the cylinder, 

*N,in AMNT>-\/Tun , aNN 



E 



n,k z 



(B.7) 



(B. 



We next derive the specific form of the T matrix co- 
efficients by matching the boundary conditions for the 
medium, i.e. the continuity of E$, E z , and H z across 
the cylindrical surface. Plugging the explicit form of 
cylindrical harmonics ( |A.1[ ) and (B.3) into these condi- 
tions we obtain two sets of four linear equations for the 
expansion coefficients. Using B = — i(c/cj)V x E and 
H = B//i we can write the system of equations for reflec- 
tion and transmission amplitudes in case of the incident 
magnetic waves in the form, 



( A MM \ 

rpMM 

\ rpNM I 

V / n,k z 



( ±qJ n (qR) \ 
J' n {qR) 




V 



(B.9) 



^UqR) J 



with the matrix 



J' n {qMR) 



n k . 



JniquR) 



(qR) 



-H' n {1) (qR) 



nk~c f / 
qRu n 



(1) 



(qR) 



7ik z 





^j^J n (qNR) 
^amj n (q N R) 

SrJ'niqNR) 



nk.c ttV 
' qRu n « 



\ 

W (qR) 



-%Hg\qR) 
-H'^\qR) 



(B.10) 



For the incident electric waves the linear equations are 



(B.ll) 



/ A MN ^ 




( o \ 


rpMN 




imuqR) 






tiMqR) 


rpNN 


n,k z 


K J'n(qR) J 



The solutions to these sets of equations (B.9 1 and (B.ll ) 



are provided in Eqs. (30)-(37). 



Appendix C: Small R expansion of the T operator of 
the cylinder 



In order to derive Eqs. (38) and (39), we need the ex- 



pansion of the T operator in terms of uR/c. For a cylin- 
der made of isotropic material with magnetic permeabil- 
ity and dielectric permittivity e(w), we find for the 
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limit R < {6, c/w}, 



rpNN _ 



( £ -i)(^-i)M/c) 2 



(C.l) 



where k z = k 
Eqs. (f22l) and 



d23b 



Substitution of the these forms into 



yields Eqs. (38) and (39 1 



T T = -i(n-l)(k 2 z - l)M?/c) 2 , (C.2) 

4 Appendix D: Leading term of T operator for 

w/c > 7? > 5 

jiNN _ rpNN = »7T fc 2 (^+ l)(g- 1) + Qi -!)(£ + !) 

~ 1,fe * 4 (e + 1)(m + 1) For w / c > # > <*) tne leading term of T operator is 

x (luR/c) 2 -^M* e l ement which has then the following form, 

(C.3) 



rpMM _ rp 
1 l,k z — 1 -l,k z 



mm _ ™ ^0 - + 1) + (m + l)(e - 1) 



lim T \ N 



4 (e + l)(/i + l) , 

x ( W i?/c) 2 , ^ + 2 ^ + (1 _ fc ; )(2t 2 + ^K/c) + 24 lQ S [V 1 - k *" R W 

(C.4) (D.l) 



rpMN _ rpNM rpM n _ rpNM _ ™ {e/i-l)k z where 1e w 0.577 is the Euler-Mascheroni constant. It 

i i,fe i -i,fcj i -i,fe^ o , , i\ can be numerically shown that other elements of T matrix 



2 (e + l)(/*+l) 
x (wi?/c) 2 , 

(C.5) 



are negligible. 
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